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We report a theoretical analysis of parametric electron 
pump through a quantum dot in the Kondo regime. In the 
adiabatic regime, we have derived the expression for pumped 
current in the Kondo regime using non-equilibrium Green's 
function. The pumped current versus different system pa- 
rameters such as gate voltage, pumping amplitude, as well as 
the phase difference between two pumping forces are calcu- 
lated and interesting physics are revealed. 
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The general physics of parametsiG-electron pump has 
been the subject of recent studieslirilj It is in particular 
inspired by theprecent experiment of Switkes et al. In 
this experimentH, the pumped current through an open 
quantum dot is driven by two gates with oscillating volt- 
ages controlling the deformation of the shape of the dot. 
The pumped DC voltage Vdot is measured to vary with 
the phase difference between the two gate voltages, 
and is antisymmetric about (j) = ti. At low pumping 
amplitude the experimental data gave Vdot ~ sin (j). In 
strong pumping regime, the dependence of Vdot on ip be- 
comes non-sinusoidal showing Vdot(O) ^ 0, whereas keep- 
ing Vdot (tt) ~ for all pumping strength. Many of these 
experimental findings have been explained theoretically. 
However, apart from Ref. to date most of the the- 

oretical investigations of parametric pumping have as- 
sumed single electron approximation. It would be inter- 
esting to see how the strong electron electron interaction 
modifies the pumped current. For this purpose, we report 
in this paper a theoretical analysis of the parametric elec- 
tron pump through quantum dot in the Kondo regime us- 
ing adiabatic theory. Our results indicate, in the Kondo 
regime, that the general behavior of the pumped cur- 
rent is similar to that of the conductance. Above the 
Kondo temperature, as one scans the gate voltage Vg we 
found two peaks in the pumped current corresponding 
to the resonant tunneling peak and Coulomb charging 
peak. When the temperature is below the Kondo tem- 
perature, a new peak in the pumped current starts to 
emerge at = — [//2 in the middle of the resonant peak 
and Coulomb charging peak. At zero temperature, the 
pumped current has a broad peak at Vg = —11/2 which is 
the superposition of these three peaks. This is very dif- 
ferent from the noninteracting case where there is only 
one peak in the pumped current. In the Kondo regime, 
we found that as one varies the pumping amplitude, the 
pumped current increases quadratically for small ampli- 
tude and then scales linearly with the pumping ampli- 



tude. Our result also shows that the pumped current is 
antisymmetric about </> = tt and is a nonsinusoidal func- 
tion of <j) for large pumping amplitude. Our result sug- 
gests that the Kondo signature can also be found in the 
pumped current which can be checked experimentally. 

We consider a 2D quantum dot with leads connected to 
the dot through narrow constrictions controlled by gate 
voltages. Since the threshold of electron propagation in 
the constriction maybe lower than that in the lead, the 
constrictions act like a double barrier whose height can 
be tuned by two gate voltages. The cyclic variation of 
these two pumping gate voltages allow the parametric 
electron pumping through the quantum dot. To simply 
the calculation, we use the one dimensional double bar- 
rier potential to model the quantum dot. 

To analyze parametric quantum pumping, we make 
use of the noncquilibrium Green's function method. 
Using the distribution function, the total charge in 
the system during the pumping is given by Q{x,t) = 
-ieJ{dE/TT){G<iE,{xlt)}))^^ where G< is the lesser 
Green's function in real space, x labels the position, and 
{X{t)} describes a set of external parameters which fa- 
cilitates the pumping process. is related to the r|e^ 
tarded and advanced Green's functions C" and G'°,t£l 
G<{E, {X}) = -f{E)[G-{E, {X}) ~ G''{E,{X\)\ where 
the retarded Green's function in real space is given by 

GHE, {X\) = i (1) 

In Eq.(|l|), S*" is the self-energy and is a diagonal ma- 
trix describing the variation of the potential landscape 
due to the external pumping parameter X. In order 
for a parametric electron pump to function, we need si- 
multaneous variation of two system parameters Xi{t) = 
Xw + XipSii\{ujt) and X2{t) = X20 + X2pS\n.{ijjt -f 4>). 
Hence, in our case, the potential due to the gates can 
be written as T^; = -''^i^i + ^2^2, where is poten- 
tial profile for each gate. If the time variation of these 
parameters are slow, i.e. for X{t) = Xq + SXsmiiot), 
then the charge of the system coming from all contacts 
due to the infinitesimal change of the system parameter 
{6X 0) is 



dQ{t) = ^9x.Tr[Q(a;,0] 5X,{t) 



(2) 



where Tr[..] is over the positions. It is easily seen that 
the total charge in the system in a period is zero which 
is required for the charge conservation. To calculate the 
pumped current, we have to find the charge dQa passing 
through contact a due to the change of the system pa- 
rameters. Using the Dyson equation dxiG^ = G^AiG^, 
Eq.(0) becomes. 
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dQ{t) = ^-^YlJ dETrlCAjG'' ~ c.c]f{E)SXjit) 

- ^ I dE{dEf) J2 nG'-rG^A,]sx,{t) (3) 

j 

where we have used the fact that G'^rG" = - G°) 
and r = Ta is the hne width function. So we obtain 

dQa{t)^^ J d£;(a£;/)^Tr[G'^r„G''A,]«, (t) (4) 



Furthermore, the current flowing through contact a 
due to the variation of parameters Xi and X2, in one 
period of time, is given by 
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dt dQa/dt 



(5) 



where r = 27p/ij is the peijiod of cycHc variation. In terms 
of injectivitytJ given byll3 



dNa 
'dX. 



27r 



Eq.(||) reduces to the famiUar formula! 

'dNc, dXi dNc dX2 



la^ / dt 



dXi dt dX2 dt 



(6) 



(7) 



Note that Eq. (|^) is a general expression applicable to the 
case of interacting and noninteracting systems as long as 
the retarded Green's function is known. 

For the transport in the Kondo regime, we consider the 
following Hamiltonian H = Hq + Hj + Ht with, 



Hq = Y^ ekaC+^^Gkca + Y^[Em + Vg]dl„^d^r 



and 



(8) 



(9) 



(10) 



where C\^^ is the creation operator of lead a and 
is the creation operator of the scattering regime at en- 
ergy level m. We have applied the gate voltage Vg to 
control the energy level in the scattering region. For this 
Hamiltonian Fq, defined in Eq.(|4|) is given by {Ta)mn — 
2tt J2k Tkam'^kanS{E - Eka) ■ There are many , 



iaches 

to treat scattering problem in Kondo regimeEj~0. We 
find it is convenience to use-the perturbation scheme pp^- 
posed by Levy Yeyati et alEJ and Kajueter and Kotliaro. 
In this approach, the retarded Green's function is given 
by 



G'' = 



E- H 



(11) 



''lead 



where is the self-energy due to the coupling between 
the scattering region and leads. The effect of strongly 
elej?t||njj-(j;lectron interaction is included in the self-energy 



K{E) = Un 



1 - B^liE) 



(12) 



where Sg is the self-energy due to the second order con- 
tribution in U , 



tip 
8^ 



dEidE2dE3 



E + E3-Ei-E2+iS 



[G>{E,)G>iE2)G<iE3) - G<{E^)G<iE2)G>{E3)] (13) 



where G^ = 1/{E- Ho ~^,,,) and G< = -/(Gg-Gg). 
Here for simplicity, we have only considered a partic- 
ular energy level Eq and used the wideband limitE^I. 
The coefficients A and B in Eq.(|l2|) are determined by 
the solutions in two limiting cases: large energy limit 
and atomic limited, from which we have A = [n{l — 
n)]/[no(l - "0)] and B = [{1 - 2n)]/[rio(l - no)U] with 
n = — J dE f {E)lmG^ / TT is the physical particle number 
and no — ~ j dEf{E)lmGo/iT is the fictitious particle 
number. This scheme gives a good description for the 
case of half filling. Away from that, one must replace Ho 
in Gg and G^ by a self-cpoaistent Hamiltonian Hf,f / and 
use the Friedel sum ruleEjCj 



n= arctan 

2 TT ^ IniS'- ^ 



(14) 



The self-consistent solution of Eqs.(|ll|), (^, and ( |14| ) 
determines the self-energy YZ. which will be used in the 
calculation of pumped current. We now apply Eq.(j^) 
to calculate the pumped current in the Kondo regime. 
The double barrier structure is modeled by potential 
U{x) = Xi5{x + a/2) + X25{x - a/2) where Xi and 
X2 are barrier heights which varys in a cyclic fashion 
to allow the charge pumping. In particular, we set 
Xi = vo -\- VpS\Q-{Ljjt + (pi) with (f>i = and 02 = <!>■ 
We will fix the units by setting h = 2m = 1 in the fol- 
lowing analysis. For the GaAs system with a = lOOOA, 
the energy uint is i? = 56/16^. We will also fix the on 
site potential U = b which is much smaller than the 
level spacing in the quantum dot, frequency uj = 1, the 
barrier height vo = 79.2, and phase difference = 7r/2 
(unless specified otherwise). Finally, the energy of in- 
coming electron is chosen to be in line with a resonant 
level Eo when Vg = 0. In Fig.l we present the transmis- 
sion coefficient versus gate voltage (which controls the 
levels in the quantum dot) at different temperatures T . 
When the temperature is higher than the Kondo temper- 
ature Tk = 0.02 (dashed line in Fig.l), we see two peaks: 
resonant tunneling peak at Vg = (for E = Eo) and 
the Coulomb charging peak at Vg = —U {E = Eo + U). 
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At low temperatures below T^, the co-tunneling process 
leads to a new peak, the Kondo peak, at the Fermi 
level. As the temperature is lowered, the peak height 
of Kondo peak increases and the dip between resonant 
peak and charging peak diminishes. At zero tempera- 
ture, the broad peak at Vg = ~U /2 in Fig.l is the super- 
position of these three peaks. For Vg > Q or Vg < —U, 
the transmission coefficient is almost temperature inde- 
pendent. Note that the peak heights (near Vg ^ and 
Vg = —U) are asymmetric about the Vg = —U/2. This 
is because the linewidth function F depends on energy 
or in our case depends on the gate voltage. Fig. 2 de- 
picts the pumped current as a function of gate voltage at 
different temperatures and for different pumping ampli- 
tudes. Generally speaking, the pumped current follows 
similar pattern of the transmission coefficient at different 
temperatures due to the fact that the pumped current is 
proportional to the density of states of the system which 
also manifests in the transmission coefficient. We see 
that as the pumping amplitude becomes larger, the ratio 
/p(T = 0,vg = -[//2)//p(T = l.^Tk.Vg = -U/2) be- 
comes smaller (see Fig. 3 for further discussion); at T ^ 
two resonant peaks become broader and move away from 
each other; the Kondo peak at T = becomes broader 
and flattened. We also notice that the pumped current 
increases as the pumping amplitude increases. In par- 
ticular, as the pumping amplitude increases, the peak 
height of resonant states at £^o and E'o + C/ increase much 
faster than that of the Kondo peak and when Vp — O.Iwq 
they have almost the same height. In Fig. 3, we plot the 
pumped current versus relative pumping amplitude Up/uo 
at two different gate voltage: one at Vg = —U/2 and the 
other near the resonant level when Vg = —0.5. At the 
Vg = —U/2, the dependence on the relative pumping am- 
plitude shows the expected quadratic behavior for small 
amplitude since the pumped current is biliacar in pump- 
ing amplitude in the weak pumping regimcfi. For larger 
amplitude Vp/vQ > 0.03 it is almost linear with different 
slopes depending on temperatures. The slope is smaller 
at higher temperature. For the gate voltage near the 
resonant level, the pumped current has similar behav- 
ior except that it is not sensitive to the change of tem- 
perature. Fig. 4 displays the pumped current as a func- 
tion of phase difference between two pumping forces 
for different pumping amplitudes. The pumped current 
is antisymmetric about the phase difference 4> = n. In 
the weak pumping regime (vp = O.OIuq), the pumped 
current shows the sinusoidal behavior and peaked at 
(j) = 7r/2. This is because in the weak pumping regime, 
the pumped current is bilimcar in the pumping amplitude 
and proportional to sin (/).a In the strong pumping regime 
(vp = 0.1 Wo), we start to see nonsinusoidal behavior as 
higher order terms of pumping amplitude come into play. 
The maximum pumped current occurs approximately at 
(j> = O.Gtt (see Fig. 4a and Fig.4b)_ Similar nonlinear be- 
havior is also seen experimcntallyd although the physical 
origin may be different. 

In summary, we have studied the parametric electron 



pumping through a quantum dot in the Kondo regime us- 
ing a nonequilibrium Green's function theory. We found 
that the behavior of the pumped current is closely re- 
lated to the conductance. As one varies the pumping am- 
plitude, the pumped current increases quadratically for 
small amplitude and then scales linearly with the pump- 
ing amplitude. Because of the resonant nature of the 
pumping, the pumped current shows nonsinusoidal de- 
pendence on the phase difference of the pumping param- 
eters. In this paper, we have used the adiabatic theory 
to calculate the pumped current. This theory is valid 
in the low frequency regime and can not account for 
the anomaly at (f) — tt found experimentallytl. At finite 
frequency, one must usCj-the real space nonequilibrium 
Green's function mcthodEil to calculate the pumped cur- 
rent. 
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FIG. 1. The transmission coefficient versus gate voltage at 
different temperatures. 

FIG. 2. The pumped current versus gate voltage for dif- 
ferent pumping amplitudes Vp. Main figure: Vp = O.lvo; left 
inset: Vp = O.Olwo; right inset: Vp = O.OSiiq- 

FIG. 3. The pumped current versus relative pumping am- 
plitude at different gate voltages. Main figure: Vg — —2.5; 
inset: Vg = —0.48. 

FIG. 4. The pumped current versus phase difference at 
different temperatures, (a). Main figure: Vg — —2.5 and 
Vp = O.Olwo; inset: Vg — —2.5 and Vp = O.Ivq. (b). Main 
figure: Vg = —0.48 and Vp = O.Olwo; inset: Vg = —0.48 and 
Vp — O.lvo- 
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